Lattice QCD study of the .s-wave nn scattering lengths in the 1 = and 2 channels 
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The s-wave pion-pion (nir) scattering lengths are computed below the inelastic threshold by the 
Liischer technique with pion masses ranging from 240 MeV to 463 MeV. In the Asqtad-improved 
staggered fermion formulation, we calculate the tvk four-point functions for the 7 = and 2 channels 
with "moving" wall sources without gauge fixing, and analyze them at the next-to-leading order in 
the continuum three-flavor chiral perturbation theory. At the physical pion mass, we secure the s- 
wave 7T7T scattering lengths as m^a 1 ^ = 0.214(4) (7) and m w al= 2 = -0.04430(25) (40) for the 7 = 
and 2 channels, respectively, where the first uncertainties are statistical and second ones are our 
estimates of several systematic effects. Our lattice results for the s-wave -kit scattering lengths are 
in well accordance with available experimental reports and theoretical forecasts at low momentum. 
A basic ingredient in our study for the 7 = case is properly incorporating disconnected diagram. 
These lattice computations are carried out with the MILC 2 + 1 flavor gauge configurations at two 
lattice spacings a ~ 0.15 and 0.12 fm. 

PACS numbers: 12.38.Gc, 11.15.Ha 
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I. INTRODUCTION 

The research on the tttt scattering is a basic and classi- 
cal subject in the field of strong hadronic interactions. 
Its handleability and simplicity essentially stem from 
the pseudo Nambu-Goldstone boson nature of pion, a 
natural aftermath of the spontaneous chiral symmetry 
breaking in quantum chromodynamics (QCD), which im- 
poses rigid constraints on the 7T7t low-energy interactions. 
Moreover, the s-wave 7T7T scattering lengths vanish in the 
chiral limit when the momentum of the pions approaches 
zero. Since these quantities stand for a sensitive probe 
of the chiral symmetry breaking generated by the quark 
masses, the lattice QCD study, an objective of this paper, 
is a non-perturbative method in an effort to comprehend 
the low-energy nature of QCD. 

With small pion masses and low-momenta, the s-wave 
7T7T scattering lengths can be solely predicted at leading 
order (LO) in chiral perturbation theory (xPT) [l|. The 
next-to-leading order (NLO) and next-to-next-to-leading 
order (NNLO) corrections in the chiral expansion 0t3 
lead to pcrturbative deviations from the LO and involve 
both computable nonanalytic contributions and analytic 
terms with some unknown low-energy constants (LEC's), 
which can be obtained from lattice simulations or exper- 
imental measurements. 

A combination of some experimental and theoretical 
inputs from CGL @, H|, along with the Roy-equation 0, 
H[ , produced a precise result of the s-wave 7T7t scattering 
lengths. Zhou ct al. studied the pole structure of the 
low-energy 7T7r scattering amplitudes using a proper chi- 
ral unitarization method in addition to the crossing sym- 
metry and low-energy phase shift data, and estimated 
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the s-wave 7T7t scattering lengths [7|. K. Sasaki and N. 
Ishizuka found that the scattering phase can be obtained 
directly from the 7T7t wave function Q. Guo et al. pro- 
vided a reliable and solid estimation of all part of the 
C(p 6 ) calculation 0, and some resonance contributions 
were added in to the former phenomenological calcula- 
tions 0, 3 , and obtained the slight differences with re- 
spect to previous results in Refs. [!, H[. Using the NLO 
SU(2) unitary chiral perturbation theory to examine the 
7T7r scattering, Albaladejo and Oiler obtained a good re- 
production of the s-wavc 7T7t scattering lengths |10j . 

In conjunction with the strict xPT constraints in the 
analysis, the considerably improved accuracy for the s- 
wave 7T7T scattering lengths has been obtained from the 
experimental measurement of the semileptonic K e 4 de- 
cay by E865 [ll|. With the independent experimental 
uncertainties and different theoretical inputs [H, [3, [Hj], 
the NA48/2 high-precision analyses of the K e 4 and 
decays [13l - ll6l | gave rise to the complementary informa- 
tion on the s-wave 7T7t scattering lengths [l5| . All of these 
theoretical (or phenomenological) predictions and exper- 
imental determinations are consistent with each other. 

Lattice studies on the 7T7t scattering have been con- 
ducted in quenched QCD by various groups [l7l - [23| . The 
full lattice study of the s-wave I = 2 irir scattering length 
was first carried out by CP-PACS [24|. Fully-dynamical 
computation of the 1 = 2 irir scattering was explored 
by NPLQCD with the domain-wall valence quarks on a 
fourth-rooted staggered sea quarks [25|, [2(| . Using the 
Nf = 2 maximally twisted mass fermion ensembles, Xu 
et al. employed the lightest pion mass at that time, con- 
ducted an explicit check for lattice artifacts [13] • With an 
anisotropic Nf = 2 + 1 clover fermion discretization, the 
I = 2 7T7T scattering phase shift is calculated by NPLQCD 
to determine all the threshold parameters [2g. More- 
over, efforts were made to first secure the <i-wave 1 = 2 
7T7T phase shift in some nice works by HSC [2j| [3(| . Us- 
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ing overlap fermion formulation, Yagi et al. examined 
the consistency of the lattice data with the NNLO xPT 
prediction after correcting finite volume effects [31| . 

Nevertheless, only a couple of lattice studies in the 
1 = channel are reported so far, whose computations 
are hindered by the so-called "disconnected diagram". 
Using the quenched approximation, Kuramashi et al. 
carried out the pioneering work for isospin-0, however, 
the vacuum diagram was disregarded assuming that vac- 
uum amplitude remains small for large t [l8|. Addition- 
ally, for the rectangular and vacuum diagrams, quark 
loops are required to make the scattering amplitudes uni- 
tary, otherwise, the basic part of the physics is lost due 
to quenched approximation [l7| . Liu conducted the first 
full QCD calculation for the 1 = channel including 
the vacuum graph, however the error of the extracted 
scattering length is remarkably large due to the usage 
of big pion masses (small one is 430 MeV) [32|]. With 
the presence of the vacuum diagram, we have attempted 
to crudely calculate the tttt scattering for isospin-0, and 
made a first lattice calculation for l^°(fi), which is a 
LEC appearing in the xPT expression of the tttt scat- 
tering length for isospin-0 [33j . Nonetheless, we used 
the partially quenched QCD to save computational cost, 
and worked with large quark masses (33|. Moreover, the 
statistical errors are underestimated since we only con- 
sidered the primary one |33| . Furthermore, we neglected 
the obvious oscillating term due to the staggered scheme. 
We understood that the statistical errors for the ratio of 
vacuum amplitude grow as e 2rriwt [H). Consequently, us- 
ing the small quark mass is very important for the 1 = 
channel. As presented later, our lattice results will indeed 
quantitatively confirm this argument, and we acquire the 
good signals of vacuum diagram for the lattice ensembles 
with small pion masses. 

To overcome the Maiani- Testa theorem (3f|, people 
usually calculate the energy levels of two-(many-)particlc 
system enclosed in a torus, and its scattering amplitudes 
can be recovered [36l - l49| . In this work, Liischer's tech- 
nique [36U38j is employed to extract the scattering phase 
shift with the lattice-calculated energy eigenstates. 

We here use the MILC gauge configurations [50|, [H| 
with the 2 + 1 flavors of the Asqtad-improved staggered 
dynamical quarks [52( to compute the s-wave 7T7T scatter- 
ing lengths. The technique of the "moving" wall source 
without gauge fixing [53| first introduced in Refs. pUGil 
are exploited to calculate all the four diagrams classified 
in Refs. [IlHlH, and special effort is payed to the discon- 
nect diagram. Our lightest pion mass is about 240 MeV, 
which is lighter than those of the former lattice studies 
on the tttt scattering and enables us to further explore 
the chiral limit. Consequently, the signals of vacuum di- 
agram are remarkably improved. Moreover, due to the 
nature of staggered fermion, our computations are au- 
tomatically precise to 0(a 2 ) 17 1. Additionally, we used 
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-0.04430(25)(40); I 
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3.27(.77)(1.12), 



where a 1 ^ 2 denote the s-wave tttt scattering lengths in 
the 1 = 2 channel and l^ 2 (fi) is a LEC evaluated at 
the physical pion decay constant. These results are in 
well agreement with the experimental measurements and 
theoretical (or phenomenological) determinations as well 
as previous lattice calculations. Most of all, we obtain 

mW=* = 0.214(4)(7); l T =° = 43.2(3.5)(5.6), 

which are in fair accordance with the experimental re- 
ports and theoretical (or phenomenological) predictions, 
and significantly improve our former study in Ref. [33j . 

The paper is organized as follows. In Sec. |TT] we will 
review the basic formalism for the calculation of the s- 
wave tttt scattering [37l |38| . The simulation parameters 
and our concrete lattice calculations are shown in Sec. lIIII 
We will give the results of the lattice simulation data in 
Sec. IIV1 fitting analyses in Sec. [V] and chiral extrapo- 
lation along with the comparisons of different results in 
Sec. IVI1 Finally, a summary of our conclusions and out- 
looks are arrived at in Sec. IVIII The compact continuum 
three-flavor %PT forms at NNLO for the tttt scattering 
lengths are courteously dedicated in Appendix 



II. METHOD 

On the basis of the original derivations and notations 
m Refs. [TM3, we reviewed the indispensable formulas 
for the lattice QCD evaluation of the s-wave tttt scat- 
tering lengths in a torus. The formulas and the no- 
tations adopted here are actually the same as those in 
Refs. [33l l54j|. But, to make this paper self-supporting, 
all the essential parts will be reiterated subsequently. 

Let us review scattering of two Nambu-Goldstone pi- 
ons in the Asqtad-improved staggered fermion formalism. 
For the s-wave tttt scattering, only the isospin 1 = and 
2 channels are permitted owing to Bose symmetry. We 
build these tttt isospin eigenchannels using the following 
interpolating operators 0, [l9[ 



ot=°(t) = 



J={.-(^+(t + l 
7r°(i)7r°(i + l)}, 

TT + (t)TT+(t+l), 



) + TT + (t)TT-(t+l)- 



(1) 



the continuum three-flavor xPT at NLO to extrapolate 
our lattice-measured tttt scattering lengths to the physi- 



with the interpolating pion operators denoted by 
7T + (t) = J(x,t)7 5 u(x,t), 

X 

7T~(t) = ^u(x,t)7 5 d(x,t), 

X 

= ^5^[«(x,t)75«(x,*)-d(x,t)7 5 d(x J t)] > 
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then we express the itit four-point function in the zero 
momentum state as 

Xi X 2 X 3 X 4 

where, to prevent the intricate color Fierz rearrangement 
of the quark lines [M, GIIj 1 we familiarly select t\ = 0, 
t-2 = 1, h = t, and U = t + 1 [Mill- 
In the isospin limit, only four diagrams contribute to 
7T7T scattering amplitudes, in Fig.Q] we show these quark- 
line diagrams, which arc identified as direct (D), crossed 
(G), rectangular (R), and vacuum (V) diagrams, respec- 
tively [M, [l9[ . The reliable evaluation of the rectangular 
diagram is challenging and the rigorous computation of 
the vacuum diagram is pretty hard [3 . 
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(a)D 



(b) C 



(c)R 



(d) V 



FIG. 1: Quark- line diagrams contributing to nn four-point 
functions. Short bars indicate the wall sources. The wall sinks 
for local pion operators are represented by open circles. 

In our former works [H, HH , we calculated these four 
diagrams via evaluating T quark propagators [IH [l9| 



(x,t)> 



< t < T- 1. 



The combination of Gf(n) which we apply for the 7T7t 
four-point functions is schematically illustrated in Fig. [TJ 
and these diagrams can be described by means of G, 



C D (ti, t,2, ^3, ti) 
C (ti, <2, ^37 ti) 
C R (ti,t2, ^3 > £4) 



^^(Tr[Gl 1 (x3,t 3 )G tl (x3,i 3 )] 

x 3 x 4 

xTr[C7 t t 2 (x 4 ,t4)Gt 2 (x4,t 4 )]>, 
^^(Tr[4(x3 J t3)G t2 (x3 J t 3 ) 

x 3 x 4 

x G t t 2 (x 4 ,t 4 )G tl (x4,t4)]>, 
^(Tr[Gj 1 (x 2 ,i 2 )G t4 (x 2 ,i 2 ) 

xG t t 4 (x 3 ,t 3 )G tl (x3,i3)]), 



1 Ficrz contributions force us to overcome the obstacle due to the 
staggered-flavor symmetry breaking [l7ll . The same problem is 
also met for the nK scattering, which is addressed by Lang et al. 
in Ref. [55| . In principle, they can be disentangled by the way 
discussed in Ref. [l7ll . but strenuous in practice. Fortunately, it 
can be trivially handled by the way introduced in Refs. (lEjUlSj. 
i.e., n meson operators are separated by a unit time slice. 



C V (ti, £ 2 , £3, ti) 



^^|(Tr[Gl i (x 2 ,t 2 )G tl (x 2 ,t 2 )] 

x 2 x 3 ^ 

xTr[G| 4 (x 3 ,t3)G t4 (x 3j t3)]) 
-(Tr[G t t i (x 2 ,t 2 )G tl (x 2 ,t 2 )) 

x(Tr[G t t 4 (x3,t 3 )G t4 (x 3 ,t 3 )]>), (2) 



where the indicated traces are conducted over color and 
the 7 5 factors are neatly removed using the Hermiticity 
attributes of the propagator G, and a vacuum deduction 
is a natural companion to the vacuum diagram [Hsj . 

The rectangular and vacuum diagrams inevitably cre- 
ate the gauge- variant noise [M, [l9| , which are neatly di- 
minished by executing the gauge field average without 
gauge fixing as we practiced in Refs. [H, [H, [EH [57l - [60| . 
In the isospin limit, the 7T7t four-point functions for the 
1 = and 2 channels can be expressed on the strength 
of four diagrams fl7l - fl9| , 



cl= (t)^(ol=°(t)\oi=°(o)y- 

G^= 2 W^(^= 2 (t)|Oi= 2 (0))= 



D 
D 



N f C, 



3NfR 



2^ 
(3) 



:tify 

s 



where the staggered-flavor factor Nf is inserted to rectif 
for the extra factor Nt in the valence fermion loops 
The four-fold degeneracy of the staggered sea quarks is 
removed by conducting the quadruple root of the fermion 
determinant [13, [M|- The fourth-root recipeis assumed to 
be able to restore the right continuum limit of QCD [6l| , 
and our results rest on this hypothesis. See Ref. [62| for 
more discussions about the fourth-root trick. 

It is customary to make use of the effective range ex- 
pansion for parameterizing the low-momentum behavior 
of the s-wave 7T7T scattering phase 5q 7 



1 1 



fccot S (k) = - + -rk 2 +0{k 4 ) 



(4) 



where a is the s-wave 7T7T scattering length, r is the effec- 
tive range, and k is the magnitude of the center-of-mass 
scattering momentum related to the energy of the 
7T7T system with total isospin / in a torus of size L by 



k 2 = \ (EL) 



k = T q > 



(5) 



here the dimensionless momentum gel. The s-wave 7T7t 
scattering length in the continuum limit is denoted by 

tan<5 (fc) 
a = hm , 

fe^o k 

which is purely elastic below inelastic thresholds. 2 Wc 
should keep in mind that the truncation of the effective 



We are only interested in the elastic region: 2m T < E^ n < 4m„-, 
where there is no Aw channel, and not up to the opening of the 
KK channel at around 1 GeV yet [63 |. where the KK channel 
contributes remarkably to the isoscalar tttt interactions I(i4l . 
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range r in Eq. (j4]) is considered as an important source of 
systematic error, which appears as Q(l/L 6 ). The 5o(k) 
can be computed by the Liischer formula [13, HH , 

kcot6o(k) = ^7r-^ 2 Z 00 (l,q 2 ), (6) 

where the dimensionless momentum q = kL / (2ir) and 
the zeta function Zqo(1; q 2 ) is formally expressed by 



2 00 (i;<7 2 ) = -j= Y] 2 1 2 ' 



(7) 



We generally compute the zeta function by the method 
discussed in Ref . [24j . Recently, an equivalent formula is 
established (65|. It allows us to avoid the subthreshold 
singularities inherent to Liischer technique [65j . 

The energy can be secured from the irir four-point 
function which manifests as 1661 



CUt) = Z^cosh[Ei„(t-\T)] 
+ (-l) t Z^cosh [Ei'„(t 



(8) 



for a large t to reduce the excited states, the terms alter- 
nating in sign is a representative feature of a staggered 
scheme (66[, and the ellipsis indicates the contributions 
from the excited states which are suppressed exponen- 
tially. In practice, the pollution due to the "wraparound" 
effects [U [I?], H?} should be taken into account. 

It should be worth to stress that, even if we project 
onto the Goldstone pions at source and sink timeslices, 
pions with all 16 stag gered- flavors still emerge at the in- 
termediate times [131 ■ However, in large t, the contribu- 
tions of non-Goldstone pions in the intermediate states 
is exponentially reduced due to their heavier masses in 
contrast with those of the Goldstone pions 

In practice, for the sake of a more intuitive presentation 
of our results, we compute the ratios 3 



R x (t) = 



C* (0,l,t,t + l) 
C,(0,^(l,t + 1)' 



X = D,C,R,a,ndV, 



(9) 

where CV(0, t) and C„(l, t + 1) are pion correlators with 
zero momentum. With the consideration of Eq. ([3]), we 
can depict the tttt scattering amplitudes which project 
out the I = and 2 isospin eigenstates as 

i?/=o(i) = R D (t) + lN f R c (t)-3N f R R (t) + ^R v (t), 



Ri=2{t) = R D {t) - N f R c (t). 



(10) 



In this work, we will employ two approaches to cal- 
culate the pion mass m^. The first method is to use 
both the point-source and point-sink operator. Never- 
theless, the point operator has big overlap with excited 
states [68[ , and in practice it is customary to use the wall- 
source operator which efficiently reduces these overlaps, 
along with a point-sink f68j . In addition we need both 
propagators to calculate pion decay constant [50|, HH . 



III. LATTICE CALCULATION 

We used the MILC gauge configurations [5(| [5l| 
with the 2 + 1 flavors of asqtad-improved staggered sea 
quarks [52j and a Symanzik- improved gluon action [69| . 
See detailed simulation parameters in Ref. [70|]. It is 
worth mentioning that the MILC gauge configurations 
are generated using the staggered formulation of lattice 
fermions [7lj with the fourth-root of fermion determinant 
which are hypercubic-smeared (HYP-smeared) [72||. As 
shown in Ref. 



73j . the chiral symmetry arc significantly 



enhanced via the HYP-smeared gauge link. 

The lattice simulation parameters of the MILC gauge 
configurations used here are epitomized in Table HI The 
simulated bare masses of light and strange sea quarks are 
denoted by ami and am s , respectively. The masses of the 
u and d quarks are degenerate, which are small enough, 
such that the physical up- and down-quark masses can be 
attained by the chiral extrapolation. The lattice spacing 
a for first three lattice ensembles is about 0.12 fm, and 
that of last three lattice ensembles is around 0.15 fm. By 
MILC convention, the lattice ensembles are referred to 
as "coarse" for a ~ 0.12 fm, and "medium-coarse" for 
a « 0.15 fm. For easy notation, it is convenient to use 
(ami,am s ) to mark lattice ensembles, e .g., "the (0.01, 
0.05) ensemble". The tadpole factors uq 17411 a re utilized 
to enhance the gauge configuration actionJHTJ, [sT| . 

To compute the irir four-point functions C^it), the 
standard conjugate gradient technique 4 is used to ac- 
quire the required matrix element of inverse fermion ma- 
trix. The periodic boundary condition is applied to the 
three spatial directions and temporal direction. We com- 
pute C W7r (t) on all the possible time slices, and collect 
them at the end of the measurement, namely, 

1 T_1 



t s =o 



After averaging the correlators over all the T possible val- 
ues of common time shift t s , as illustrated later, we found 
that the statistics are indeed significantly improved. 



3 In principle, when t <^ T/2, even if placing the periodic bound- 
ary condition in the temporal direction, the energy shift of the tttt 
system can be still roughly evaluated from these ratios. In this 
work we do not use these ratios to quantitatively calculate any 
physical quantities, nonetheless, these ratios will indeed help us 
comprehend qualitatively or intuitively some physical quantities. 



4 The conjugate gradient residual used in this work is 1.0 X 10 — 5 , 
which is smaller than that used in generating gauge configura- 
tions l50l |51| . Moreover, to avoid the potential roundoff errors 
as much as possible, all the numerical calculations are calculated 
in double precision. 
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TABLE I: The parameters of MILC gauge configurations used in the present work. The lattice dimensions are expressed in 
lattice units in the second block with spatial (L) and temporal (T) size. The gauge coupling /3 = 10/ g 2 is shown in Column three. 
The fourth and fifth blocks give the bare masses of the light and strange quark masses in terms of ami and am 3 , respectively. 
Tadpole-improvement factor uo is listed in Column six. The ratio ri/a is provided in Column seven (see Ref. [75^ for the MILC 
definition of n). The inverse lattice spacing a -1 is recapitulated in Column eight (for the (0.00484, 0.0484) ensemble, we obtain 
the value of ri/a from Ref. [7(j, then calculate a -1 ). In the last column the number of gauge configurations is given. 



Ensemble 


LxT 


P 


ami 


am s 


Mo 


n/a 


a _1 GeV 


N cf 










a « 0.12 fm 












2464f21b676m005m050 


24 s 


x 64 


6.76 


0.005 


0.050 


0.8678 


2.647(3) 


1.679l^ 


156 


2064f21b676m007m050 


20 3 


x 64 


6.76 


0.007 


0.050 


0.8678 


2.635(3) 


1.6721^ 


200 


2064f21b676m010m050 


20 3 


x 64 


6.76 


0.010 


0.050 


0.8677 


2.619(3) 


1.663±te 


200 










a « 0.15 fm 












2048f21b6566m00484m0484 


20 3 


x 48 


6.566 


0.00484 


0.0484 


0.8602 


2.162(5) 


1.373±?4 


560 


1648f21b6572m0097m0484 


16 3 


x 48 


6.572 


0.0097 


0.0484 


0.8604 


2.140(4) 


1.358±?| 


250 


1648f21b6586m0194m0484 


16 3 


x 48 


6.586 


0.0194 


0.0484 


0.8609 


2.129(3) 


1.352±?| 


200 



For each time slice, three fermion matrix inversions are 
needed corresponding to the 3 color choices for the pion 
source, and each inversion takes about 1000 iterations 
(about 2000 for the (0.00484,0.0484) and (0.005,0.05) 
ensembles) during the conjugate gradient calculation. 
Thus, totally we carry out 3T inversions on a single gauge 
configuration. As shown later, this rather big number of 
the inversions offers the substantial statistics needed to 
get the 7T7T scattering amplitudes with high accuracy. 

In practice, we calculate the pion correlators, 



T-l 



t„=o 

Cf P (*) = ^J2(0\^(t + t s )W n (t s )\0), (11) 

t 3 =0 

where it is the pion point-source operator and W„ is the 
pion wall-source operator [5(J UK ■ To simplify the nota- 
tion in this section, the summation over the lattice space 
point in sink is not written out. In this work, we will 
adopt the shorthand notation: "PP" for the point-source 
point-sink propagators, and "WP" for the wall-source 
point-sink propagators [5l| . We should stress that the 
summations are also taken over all the time slice for the 
pion propagators, and we found that the statistics are 
indeed significantly improved. This is very important to 
obtain pion mass with high accuracy. 

Overlooking the excited state contributions, the pion 
mass TOtt can be secured at large t with a single exponen- 
tial fit ansatz [HI, [zj III] 



iPP 



CT(*) = K 



,-m^T-t) 



(12) 
(13) 



where T is the temporal extent of the lattice, A^ F and 
are overlapping amplitudes. We will use these val- 
ues to estimate the wrap-around contributions [2ll . l27l . l67j 
and calculate the pion decay constant [68| as well. 



Wc should remark at this point that, in the calculation 
of the 7T7r four-point functions for the 1 = channel, 
we try our best effort to compute the vacuum diagram, 
since the other three diagrams can be relatively easily 
calculated. We found that the vacuum diagram plays a 
critical role in this correlator. 5 



IV. LATTICE SIMULATION RESULTS 

A. Pion mass and pion decay constant 

In practice, the it propagators were fit by varying min- 
imum fitting distances D m ; n , and with maximum dis- 
tance D max either at T/2 or where the fractional sta- 
tistical errors surpassed about 20% for two sequential 
time slices [50|, HH, HJ • 6 In this work, pion masses were 
secured from the "effective mass" plots for each of the 
MILC lattice ensembles, and they were strenuously opted 
by looking for a combination of a "plateau" in the mass 
as a function of D m ; n , good fit quality (i.e., x 2 /dof < 1) , 
and D m j n large enough to reduce the excited states [25l — 
[27| . The WP propagators were fit to Eq. (fT3"|) using a 
minimum time distance of 14a for the "medium-coarse" 
lattices and 20a for the "coarse" lattices, and the full co- 
variance matrix is used to compute statistical errors. At 
these large distances, the pollution from excited states is 
at most comparable to the statistical errors [68| . For ex- 
ample, Fig. [2] exhibits the results for pion masses and 
amplitudes as a function of D m i n for the (0.007, 0.05) 



5 In our previous work 1771 . we presented the detailed procedure 
to calculate the disconnected diagram for the /o(600) meson. It 
helps us a lot to implement the evaluation of the vacuum diagram 
here, especially for how to conduct a vacuum subtraction. 

6 Since the lattice data points at the largest distances contain 
relatively little information, the exact selection of large distance 
cutoff D max is not very critical |50| . [5ll , l6l| . 
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ensemble. Since the major objective of this work is to 
present the work for isospin-0, as explained later, for the 
7T7T scattering in the I = channel at this work, the sys- 
tematic error of the energy of the tttt system is pretty 
large, we can temporarily neglect the systematic effect 
for pion mass due to excited states. All of these fitted 
values of pion masses are listed in Table [TTl 



0.200 



(0.007,0.05) 



<>- 



.<><><><X> 



— 



.<> 



<> 



CD 



©CDCDCO©(DCDCD©©CDa)©CDcD©CDCDCDCD®CD 



0.198 

( 

0.196 
0.194 
0.192 
0.190 
0.188 
0.1 B6 

4 6 8 10 12 14 16 18 20 22 24 26 28 30 
D min 

FIG. 2: (color online). Pion masses (magenta octagons) and 
WP amplitudes A^ F (red diamonds) as a function of the min- 
imum time distance in the fit for the (0.007, 0.05) ensemble. 
The amplitudes have been divided by 2060. 

In our previous work |33| . we used the method de- 
scribed in Ref. [68| to extract the pion decay constant 
for the (0.0097,0.0484) ensemble [33|. In the light of 
same procedures, we calculated the pion decay constants 
for other "medium coarse" ensembles. All of these fitted 
values of pion decay constants are listed in Table [TTJ 

As a consistency check, the PP correlators were reli- 
ably measured in this work. Use these correlators, we 
can secure pion masses via Eq. (|12[) which are listed in 
last block in Table HH and these pion masses are found 
to be consistent with its counterparts extracted with WP 
propagators, which are summarized in Table [TTl 

Our fitted values of pion masses and pion decay con- 
stants listed in Table [TTl are in rather good agreement 
with the same quantities which arc computed on the same 
lattice ensembles by MILC collaboration in Refs. [70l.[76j. 
For the "coarse" ensembles, the MILC's determinations 
on pion decay constants are directly quoted [zO, |7(| . 
which are also summarized in Table [TTJ 



(X = D,C,R and V) arc illustrated as the functions of 
t for each lattice ensemble. 

The ratio values of the direct amplitude R D are quite 
close to oneness, indicating a pretty weak interaction in 
this channel. The crossed amplitude, in another aspect, 
increases linearly up to t ~ 18 for the "medium coarse" 
ensemble and t ~ 24 for the "coarse" ensemble, implying 
a repulsive interaction between two pions in the 1 = 2 
channel. In contrast, after a beginning increase up until 
t ~ 4, the rectangular amplitude demonstrates a roughly 
linear decrease up until £ ~ 18 for the "medium coarse" 
ensemble and t ~ 24 for the "coarse" ensemble, suggest- 
ing an attractive force among two pions for the 1 = 
channel. Additionally, the magnitude of the slope is sim- 
ilar to that of the crossed amplitude but with opposite 
sign. Furthermore, we observe that the crossed and rect- 
angular amplitudes have the same values at t = and 
close ones for small t. These characteristics are in good 
keeping with the theoretical predictions [TtJ ■ 

It is extremely noisy in the amplitude of the vacuum 
diagram (V) for the (0.0194, 0.0484) ensemble, while we 
can see a good signal up to t = 10, and loss of signals after 
that. These characteristics are in well accordance with 
the Okubo-Zweig-Iizuka (OZI) rule and x?T in leading 
order, which expect the disappearing of the vacuum am- 
plitude [I3-G3- For the lattice ensembles with the pion 
mass changing smaller, the signals of the vacuum dia- 
gram are becoming more and more tolerating. For the 
(0.00484, 0.0484) ensemble, the signals of the vacuum di- 
agram is already tolerating. This qualitatively confirmed 
the analytical arguments in Ref. [34J which indicates that 
the error of the vacuum amplitude grows exponentially 
as e 2m,rt . The numerical calculation of the amplitude for 
the vacuum diagram stands as our principal and distinc- 
tive accomplishment of this paper. 

The systematically oscillating behavior of the rectan- 
gular amplitude in large t is evidently observed, which 
is a typical feature of the staggered formulation of lat- 
tice fermions and corresponds to the contributions from 
the intermediate states with opposite parity (6(| , and for 
the lattice ensemble with large pion mass, this oscillating 
feature become more obvious. In contrast, for that with 
small pion mass, this feature is not appreciable, and even 
not perceptible for the MILC (0.00484, 0.0484) ensemble. 
The physical meaning of this fascinatingly oscillating be- 
havior is easily understood [66j . Nevertheless, its quan- 
titative mass dependence is not clear to us, which are 
highly needed to be further investigated. 



B. Diagrams D,C,R, and V 

As practiced in our former work [33j , the rnr four-point 
functions are robustly calculated on six MILC lattice en- 
sembles listed in Table UJ using the technique of the mov- 
ing wall source without gauge fixing [l8|,[l9|,[53j ■ In Fig- El 
the individual ratios, which are denoted in Eq. ©, R x 



C. The errors of R v (t) and R R (t) 

According to the analytical arguments in Ref. [Hf , the 
error of the ratio for the vacuum amplitude increases ex- 
ponentially as e 2 " 1 "*. Therefore, it is pretty difficult to 
secure the correct signal for large t [34|. Likewise, the 
ratio for the rectangular diagram has errors, which grow 
exponentially as e" 1 "* for large t [H]]. Our lattice data 
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FIG. 3: (color online). Individual amplitude ratios R x (t) of the tyty four-point functions computed via the moving wall source 
without gauge fixing as the functions of t for six MILC lattice ensembles. Direct diagram (red diamonds) displaced by 0.8, 
vacuum diagram (magenta octagons), crossed diagram (red squares) and rectangular diagram (blue diamonds). 
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TABLE II: Summaries of the pion masses and pion decay constants. The third block shows the pion masses in lattice units 
and the fifth block give the overlapping amplitude using WP propagators. The product of m n L is presented in Column four. 
The pion decay constants in lattice units are provided in sixth block, where the star on the superscript indicate the MILC's 
determination. Column two shows the values of pion mass in MeV, where the errors are estimated from both the error on 
lattice spacing a and statistical errors in Column three. The seven block shows the dimensionless ratio m-n/ f w , where the errors 
are estimated from the ara^ and af^. The last block shows the pion masses which are measured by the point- wall point-sink 
propagators, and only used as a consistency check. 



Ensemble 


m^(MeV) 


WP 

am™ 




4 WP 






pp 


(0.00484, 0.0484) 


240(4) 


0.17503(09) 


3.5006(18) 


1000.18 ± 1.580 


0.11767(45) 


1.4874(57) 


0.17504(09) 


(0.005, 0.05) 


268(5) 


0.15970(15) 


3.8345(48) 


770.534 ± 2.577 


0.09054(33)* 


1.7639(66) 


0.15992(16) 


(0.007,0.05) 


315(6) 


0.18868(22) 


3.7736(44) 


399.094 ± 1.393 


0.09364(20)* 


2.0149(49) 


0.18871(24) 


(0.0097, 0.0484) 


334(6) 


0.24566(18) 


3.9306(29) 


395.107 ±1.151 


0.12136(29) 


2.0242(51) 


0.24587(21) 


(0.01,0.05) 


373(7) 


0.22455(27) 


4.4910(54) 


365.595 ± 2.039 


0.09805(14)* 


2.2902(42) 


0.22447(17) 


(0.0194, 0.0484) 


463(8) 


0.34279(19) 


5.4846(30) 


315.695 ± 0.865 


0.13055(48) 


2.6258(98) 


0.34279(23) 



indeed demonstrates such dependence with the expected 
slopes. 

The magnitudes of these errors are quantitatively in 
line with these analytical predictions as demonstrated 
in Fig. gj Fitting the errors 8R v {t) and 5R R {t) by a 
single exponential fit ansatz SR (t) ~ e^ vt and 5R R (t) ~ 
e MR *, respectively, for six lattice ensembles, we extract 
the fitted values of \iy and fin, which are summarized in 
Table Hill together with their fitting ranges. 



TABLE III: Summaries of the fitted values for jiv and (ir 
in lattice units. The second and third blocks show the fitted 
values of fiy, and fiR, respectively, and Column four gives the 
time range for the chosen fit. 



Ensemble 




afiR 


Range 


(0.00484, 0.0484) 


0.3392 


0.1621 


10-18 


(0.0097, 0.0484) 


0.4956 


0.2457 


8-16 


(0.0194,0.0484) 


0.6927 


0.3487 


8-16 


(0.005,0.05) 


0.3178 


0.1513 


10-20 


(0.007,0.05) 


0.3701 


0.1895 


10 - 20 


(0.01,0.05) 


0.4463 


0.2237 


10 - 20 



From Table IIII1 we note that the fitted values of /ir 
can be compared with the pion masses m r listed in Ta- 
ble [TT1 and half of the fitted values of fiy can also be 
reasonably compared with these pion masses. We here 
have numerically confirmed the Lepage's analytical ar- 
guments [HI about the 7T7t scattering. This testifies the 
practical applicability of the moving wall source without 
gauge fixing from another point of view. Thus, we can 
reasonably assume that the vacuum amplitude remains 
small for large t. In principle, we can overlook the vac- 
uum amplitude in the rest of the analysis. However, we 
will explicitly include it for the sake of completeness of 
the lattice QCD calculation. 
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FIG. 4: (color online). The errors of the amplitude ratios 
R x (t)(X = V, R) as the functions of t for each of the MILC 
lattice ensembles. Solid lines are single exponential fits, and 
the fitting ranges are listed in Table |III] (a) The errors of the 
amplitude ratios R v (t); (b) Those of R R (t) 
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D. Ri projected onto the 1 = and 2 channels 

The ratios Ri(t) projected onto the isospin 1 = 
and 2 eigenchannels for the MILC (0.00484,0.0484) and 
(0.005, 0.05) ensembles, are demonstrated in Fig. [5] A 
decrease of the ratio Ri=2 (t) indicates a repulsive inter- 
action among two pions for the 1 = 2 channel, on the 
other hand, an increase of the ratio Ri=o(t) suggests an 
attractive interaction for the 1 = channel. In the 1 = 
channel, a dip at t = 3 for the (0.00484,0.0484) and 
t = 5 for the (0.005, 0.05) can be clearly observed, and 
its physical origin is not clear to us as well [lH Ell ■ 



1 = channel in large t is not clearly observed, which is 
a typical characteristic of the Kogut-Susskind formula- 
tion of lattice fermions (6(| . Addtionally, this oscillating 
feature is hardly noticed for the 1 = 2 channel. 

In order to present the contribution from the vacuum 
term more intuitively, we employ the yellow cross points 
to indicate the ratio Ri(t) for the 7T7t four-point function 
in the 1 = channel without the presence of the discon- 
nected diagram. From Fig. [5j we can clearly notice that 
the contribution from the disconnected diagram is only 
obvious when t > 20 for the (0.00484, 0.0484) ensemble, 
and t > 16 for the (0.005, 0.05) ensemble. 
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FIG. 5: (color online). 7?/(t) (1 = and 2) for the nn four- 
point function at zero momenta calculated by the moving wall 
source without gauge fixing as the functions of t for the MILC 
lattice (a) (0.00484,0.0484) ensemble; and (b) (0.005,0.05) 
ensemble. The cross yellow points indicate the ratio Ri(t) in 
the 7 = channel where the vacuum diagram is turned off. 

Due to the rather small quark mass of two lattice en- 
sembles, the systematically oscillating behavior for the 



E. Lattice artifact 

From Fig. [3l we observe that there exists a pollu- 
tion from the "wraparound" effects [2l], [27], [13] approxi- 
mately starting as early as at t = 12 ~ 18 for the MILC 
"medium-coarse" ensembles and t = 22 ~ 25 for the 
"coarse" ensembles. As discussed in Refs. [2ll. [27l [3ll [67| . 
one of two pions can propagate T — t time steps back- 
wards due to the periodic boundary condition in the tem- 
poral direction. This operates as a constant contribu- 
tion and deforms the 7T7T four-point functions in large t 
(especiall y a round T/2), according to the discussions in 
Refs. [2lL l27l [3ll lo7| . it is roughly suppressed by 

exp (-nijl) / exp (—2m^t) 

as compared with forward propagation of the 7T7t state. 
We can opt the fit ting ranges satisfying i max <C T/2 to 
reduce this effect 1211] . However, according to the ar- 
guments in Refs. [2l[ WK loll . if pion mass is small 
enough (e.g. the (0.00484,0.0484) ensemble), the wrap- 
around pollution can not be suppressed even for t <^ T/2, 
we should include this term for the successful fit, 7 



CUt) = C + A TW cosh [E^(t- IT)] 
+(-l)*4 nr cosh [E&{t- |T 



(14) 



where C is a constant corresponding to the wrap-around 
term. This can be easily understood by evaluating the 



contribution of two fake diagrams in Fig. 2 of Ref. [31 
and the C can be expressed as 



C 



2Ah 



(15) 



For easy notation the superscript WP in A v is omitted 
in the rest of analyses. 

In this work, we accurately extract the overlapping am- 
plitudes A n and pion masses m ff corresponding to pion 



7 It turned out that a five- parameter cosh- fit of C, A^n, E^ n , A 1 ^^, 
and yields a satisfactory result with an pretty acceptable \ 2 ■ 
Moreover, the excited states will be taken into account as one of 
the important sources of systematical error in this work. 
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correlators [zU which are listed in Table |TT1 and these 
values are sufficiently precise to estimate the wraparound 
terms with Eq. (|15|) . which arc listed in Tabic llVl 



TABLE IV: Summaries of the calculated wraparound contri- 
butions from overlapping amplitude A w and pion mass am x . 
The second block shows the wraparound contributions calcu- 
lated from Eq. Q15p . where its errors are roughly estimated 
from the statistical errors of and m^. 



Ensemble 



C 



(0.00484, 0.0484) 
(0.005,0.05) 
(0.007,0.05) 
(0.0097, 0.0484) 
(0.01,0.05) 
(0.0194, 0.0484) 



451.36(2.63) 

43.23(51) 

1.815(29) 

2.362(25) 

0.1533(31) 

0.01424(15) 



We note that Rcf. 
definition as 



[30l | has recently taken the similar 



C — 2A^ 



-m„T 



(16) 



where is defined in Eq. (fT4"|). Additionally, there 
exists the similar general form which contains the two- 
particle as well as one-particle eigenvalues, and gives a 
not rigorous proof in Ref. [79| . Furthermore, when study- 
ing K — > 7T7T decay amplitudes, there is an analogous 
fitting functional form in Ref. [56j . 

In order to comprehend this wraparound effects at a 
quantitative level, we denote a quantity 



Rwc{t) — 



C 



(17) 



which is the ratio of the wraparound pollution to the 
7T7r four-point function in the 1 = 2 channel. In fact, 
we exploited the data of the wraparound contribution C 
listed in Table [TV] and C^ 2 (t) calculated from Eq. © 
to approximately evaluate these ratios. The ratios for 
six MILC lattice ensembles are illustrated in Fig. [5] All 
of these ratios make a significant contribution and are 
approximately close to 1/2 as t approaches to T 12 as 
expected from the arguments in Refs. [2l|, [I?], [U [67j . 
We can note that as the pion mass of the lattice ensem- 
ble becomes smaller, the wraparound contribution C is 
clearly observed even at small t. For example, we can 
keenly notice the wrap-around term even as early as at 
t = 10 ~ 12 satisfying t < T/2 for the (0.00484,0.0484) 
ensemble. It is, therefore, absolutely necessary for us to 
explicitly consider the wrap-around term, especially for 
lattice ensembles with small pion masses when we extract 
the energy E of the tttt system [3~0 ]. 

To get rid of this pollution, Xu et al. [27j employed a 
derivative method and denoted a modified ratio [67| . By 
ignoring terms suppressed relative to the leading contri- 
bution, the energy shift 5E can be obtained from the 
asymptotic form of the modified ratio. To identify the 
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FIG. 6: (color online). The ratios of the wraparound terms 
to the corresponding nn four-point functions for six MILC 
lattice ensembles calculated by Eq. (|17|l . All of these ratios 
make a significant contribution and are approximately close 
to 1/2 as t a ppr oaches to T/2 as expected from the arguments 
in Refs. [SEMES Ifil. 



time-separations where the ground state dominates, Yagi 
et al. used their self-defined ratios [3lj . Moreover, Dudek 
et al. H3] recently eliminated this unwanted pollution 
term by the shifted correlator. 

In principle, we can use one of the three above- 
mentioned methods to process our 7T7t scattering data 
for isospin-2. However, for those of the 1 = chan- 
nel, there is a further complication introduced by stag- 
gered fcrmions, the oscillating term is appreciable, we 
must consider the oscillating term and modify the cor- 
responding functional forms, as a consequence, it is not 
convenient to use these methods. Additionally, for the 
pion-kaon (ttK ) scattering [53], [78| , and the tttt correla- 
tors "in flight" |30j . the wrap-around term is not a con- 
stant, and these methods are not suitable. 

Nagata ct al. solved this problem by subtracting the 
wrap-around term numerically from the obtained quan- 
tities [lH , since the lattice-measured data are sufficiently 
precise to allow such subtraction. Dudek et al. elim- 
inated this term by means of the shifted correlator for 
the 7T7r scattering at rest and weighted-shifted correla- 
tor for that "in flight" [30(. We already exhibited that 
the overlapping amplitude A^ and pion mass can be 
calculated with high accuracy, so, it is natural to bor- 
row these methods to our case [3(| Iz3 • As a consistency 
check, we numerically compared these results calculated 
from the above-mentioned methods for lattice data in the 
1 = 2 channel, and found that they are well consistent 
with each other within errors. Therefore, in this work we 
only present the results from the last method, namely, 
using equation (fl"4"|) to extract the energy E of the 7T7T 
system in a conceptually clean way. 
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V. FITTING ANALYSES 

As already explained in previous sections, we will use 
Eq. ([M]) to get the energies aE of the -kit system, which 
are inserted into the Luscher formula ([5]) to obtain the 
corresponding s-wave tttt scattering lengths. Hence, ap- 
propriately extracting the energies is a central step to our 
ultimate results. A persuasive way to process our lattice 
data is the resort to the "effective energy" plot, which is 
a variant of the effective mass plot, and very similar to 
the "effective scattering length" plot [1^, [26| . 



A. 1 = 2 channel 

In practice, tttt four-point functions were fit by vary- 
ing the minimum fitting distances D m ; n , and with the 
maximum distance D max either at T/2 or where the frac- 
tional statistical errors surpassed about 20% for two se- 
quential time slices [5(1 [EH [6l| . Additionally, the fitting 
parameter C was constrained by priors to conform the 
lattice-calculated wraparound contribution C listed in 
Tables ITVl ;80j] . For each ensemble, the "effective energy" 
plots as a function of D m j n are illustrated in Fig. [7] The 
central value and statistical error at each time slice were 
evaluated by the Levenberg-Marquardt algorithm [8lj . 
To make these fits more robust, we double-check them 
with SNOBFIT, which is a soft constrained noisy opti- 
mization [s3] - From Fig. [7j we also observed that the ef- 
fective energies have larger statistical errors near t ^ T/2 
because of the wrap-around effect as it is discussed in de- 
tail in Ref. [H|. 

We should remark that the phy sical model in Eq. (TH)) 
just include the ground state [25l - t27| . In fact, we can fit 
with the inclusion of the first excited state, and the differ- 
ence between these procedures, as well as the difference 
arising from the arbitrary choice of D max , is incorporated 
in the systematic error for aE at each time slice. 

In this work, the energies aE of the mr system in the 
1 = 2 channel were secured from the "effective energy" 
plots for each of the MILC lattice ensembles, and they 
were strenuously opted by looking for a combination of 
the "plateau" in the energy as the function of D m i n , good 
fit quality [251427} , and D m i n large enough to suppress the 
excited states. We found that the effective energy of the 
tttt system for the 1 — 2 channel have relatively small 
errors within a broad minimum time distance region. 

For the same fitting range, analogously we secured the 
wraparound term C from the corresponding "effective 
wraparound constant" plots in Fig. [5] It is worth men- 
tioning that the fitted wraparound pollution C are in fair 
agreement with the calculated wraparound pollution C 
within errors. Nonetheless, it is interesting to note that 
there exist about 1% differences, and the physical indi- 
cation is not clear to us, which are highly needed to be 
further investigated in the future work. 

The fitted values of the energies aE of the tttt system 
with isospin-2, fitting range and fit quality (x 2 /dof) are 
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FIG. 7: (color online). The "effective energy" plots as the 
functions of D m m for the tttt scattering in the 1 = 2 chan- 
nel in lattice units. The "effective energy" plots have small 
errors within a broad minimum fitting distance region. The 
estimates of the systematic uncertainty due to fitting are not 
displayed in this figure. 



tabulated in Table [V] together with the fitted values of 
the wraparound contribution C. We note that the fitted 
values of C is a statistically significant constant term de 
facto for the lattice ensembles with small pion masses. 
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FIG. 8: (color online). The "effective wraparound constant" 
plots as the functions of D m i n for the -kit scattering in the 
1 = 2 channel in lattice units. 



states is a quite small fraction of total energy of the 7T7T 
system, 8 which can be estimated from the data in Ta- 
bles [III and EI This forces us to make the rigorous mea- 
surements of both the energies spectrum of the 7T7t system 
and pion masses, and even seriously account for various 
small systematic effects to resolve the rather small differ- 
ences. We have indeed extracted the energies of the 7T7T 
system and pion masses with significantly high precision 
which arc shown in Tables [XT] and IVl 

Now it is straightforward to substitute these energies 
aE into Luscher formula (J5]) and secure the relevant 
s-wave scattering lengths a^ 2 , where we plugged the 
pion masses in Column three in Table HU The center-of- 
mass scattering momentum k 2 is computed by Eq. ([5]) 
with pion masses given in Table [TTJ However, to get 
rid of the scale-setting uncertainties, it turns out to be 
more customary to adopt the dimensionlcss quantity: 
m^a 1 ^ 2 [25|, |26| | . All of these values for each lattice en- 
semble are summarized in Table [Vj where the statisti- 
cal errors of k 2 are calculated from the statistical errors 
of aE and am T , and its systematic errors are only esti- 
mated from the systematic errors of aE. Likewise, the 
statistical errors for m^a 1 ^ 2 are computed from the sys- 
tematic errors of k 2 and am^, while its systematic errors 
are estimated from the systematic errors of k 2 and the 
subsequently-mentioned two finite volume effects. 

Since the periodic boundary condition is imposed in 
the spacial directions of the lattice, there has an expo- 
nentially small finite volume (FV) correction to the s- 
wave 7T7T scattering length in the 1 = 2 channel, which 
has been determined in the vicinity of the threshold in 
Ref. [H[. The consequent finite volume correction Apy 
is provided here as [83[ , 



J=2 



)l = (m^a 1 ^ 2 )^ + A FV 



(18) 



where 



A FV = 



1 



213/2^5/2 ^ 

17 1 

x 1 



I nlm-n-L 



169 1 
8 \n\m n L ' 128 \n\ 2 mlL 2 



■ + 



+0(zr 3 )L(i9) 



here X^nez 3 indicates a summation without n = 0. Using 
this formula, we compute the finite volume corrections to 



m^a^ , which are listed in Table IVIl where we insert the 
values of m T L and m^/ f n listed in Table ITll 

From Table |VH we note that these corrections are more 
and more important for the lattice ensembles with pion 
masses smaller and smaller (83|. Since we use the lattice 
ensembles with small pion masses, we should consider 



Additionally, we clearly found that these fitted values of 
C are close to our estimated values listed in Table IIVI as 
already noticed in Ref. (30j . 

It is well-known that the interaction among two pions 
in the 1 = 2 channel is pretty weak such that the energy 
difference between the interacting and non-interacting 7T7t 



8 In this work, the ratio of the energy shift to total energy is 
about 2%. For other lattice studies I17H31H , it is actually close 
to this number. On the other hand, this ratio for the / = 
channel is around 5% [l8l . IT9I . It J , and the ratio of this work is 
approximately to 5% as well. 
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TABLE V: Summaries of lattice results of the s-wave scattering lengths for the 7 = 2 channel. The second block presents 
the energies in lattice units, where the first uncertainties are statistical, the second ones are the estimates of the systematic 
uncertainties due to fitting, Column three shows the fitted values of the wraparound term C, Column four indicates the 
time range for the chosen fit, and Column five gives the fit quality j^/Aof. Column seven gives the center-of-mass scattering 
momentum k 2 in GeV, and Column eight presents the product of the pion mass and scattering length: m^a 1 ^ 2 , where the first 
uncertainty is statistical and, the second one is systematic. 



Ensemble 


aE 


C 


Range 


X 2 /dof 


fc 2 [GeV 2 ] 




(0.00484, 0.0484) 


0.35520(25)(20) 


438.05(1.94) 


14- 


-24 


4.21/6 


0.00167(10)(7) 


-0.0915(52) (35) 


(0.005,0.05) 


0.32424(35) (33) 


42.44(29) 


16- 


-32 


10.6/12 


0.00220(21)(15) 


-0.125(H)(8) 


(0.007,0.05) 


0.38606(44) (37) 


1.776(17) 


16- 


-32 


11.5/12 


0.00444(31) (20) 


-0.167(10)(7) 


(0.0097, 0.0484) 


0.50087(41) (38) 


2.320(14) 


14 


-24 


6.1/6 


0.00437(25) (17) 


-0.167(9)(6) 


(0.01,0.05) 


0.45648(56)(41) 


0.1493(18) 


18- 


-32 


8.7/10 


0.00464(48) (26) 


-0.209(19)(10) 


(0.0194,0.0484) 


0.69392(40) (34) 


0.01414(11) 


13- 


-24 


6.6/7 


0.00527(35) (22) 


-0.277(16) (10) 



TABLE VI: Summaries of the finite volume corrections Afv- 
Column two shows the finite volume corrections to the 1 = 2 
7T7r scattering length, and Column three gives ratios of the fi- 
nite volume corrections to the corresponding statistical error. 
Here we use the pion masses, pion decay constants and m n L 
values listed in Table [TT1 



Ensemble 



Ratios 



(0.00484, 0.0484) 

(0.005,0.05) 

(0.007,0.05) 

(0.0097, 0.0484) 

(0.01,0.05) 

(0.0194,0.0484) 



0.000258 
0.000341 
0.000625 
0.000529 
0.000453 
0.000257 



0.049 
0.031 
0.060 
0.062 
0.023 
0.016 



these effects, although they are slight, and never more 
than 7% of the corresponding statistical errors [27j . 

Another important finite volume effect stems from 
effective range approximation, k cot 5(k) = 1/a^ 2 + 
\rk 2 [26|. While the dependence on the effective range 
r is small, and the range truncation actually leads to the 
correction at 0(L~ 6 ) in Liischer formula ^ [26[ . In prac- 
tice, we compute this correction for each lattice ensemble 
as suggested in Ref. [26j |. 

These two finite volume corrections have been also 
added in quadrature to the systematical errors listed in 
Table |V] Other sources of systematic uncertainty like: 
isospin violation, finite volume effect due to the fixed 
global topology, pion mass correction [26|, HH, etc. are 
believed to be very small or we currently do not have 
enough computational resources to fulfil it. These effects 
should be incorporated into the more sophisticated lat- 
tice computation at some points in the future. 



should stress that when using physical fitting model (fl"4")) 
to extract the desired energies aE of the 7T7T system, we fix 
the fitting parameters of wraparound contribution C with 
the estimated values listed in Tables IIVI 9 In practice, 
the 7T7T four-point functions were fit by altering the min- 
imum fitting distances D mm , and putting the maximum 
distance D max either at T/2 or where the fractional sta- 
tistical errors exceeded about 20% for two sequential time 
slices EH . The "effective energy" plots as the func- 
tions of D min are illustrated in Fig. [§] The central value 
and statistical error at each time slice were evaluated by 
Levenberg-Marquardt method [sT|] . To make these fits 
robust, we double-check them with SNOBFIT [12]. We 
do not show the result of (0.0194, 0.0484) ensemble in 
Fig. [9] since it is too noisy. 

For each lattice ensemble, the energies aE of the irir 
system for the 1 = channel are secured from the "ef- 
fective energy" plots, and chosen by looking for a combi- 
nation of a "plateau" in the energy as the function of 
D m in j a good confidence level and D mm large enough 
to suppress the excited states [25T - l27j ] . In addition, as 
performed for isospin-2, we approximately estimate the 
systematic errors owing to fitting [25l - l27j . which are not 
displayed in Fig. O 

The fitted values of the energies aE of the 7T7t system, 
fit range and fit quality (% 2 /dof) are summarized in Ta- 
ble |VlTl The fit quality x 2 /dof is reasonable for the 1 = 
channel. It is well-known that the disconnected term give 
rise to the considerable fluctuations to the 7T7t four-point 
function, and it is pretty hard to reliably calculate this 
term. In reality, only the lattice ensembles with small 
pion mass have a good signal. From Fig. [9] we found 
that, for the (0.097,0.0484) and (0.01,0.05) ensembles, 
the plateaus are not too obvious, so in this work, we 
don't include these results. 



B. 1 = channel 

As already performed for the 1 = 2 channel, we ana- 
lyze our lattice data with the "effective energy" plot. We 



9 For rectangular (R), and vacuum (V ) diagrams, there is no the 
wraparound pollution. So the wraparound contribution for the 
/ = channel is the same with that in the 1 = 2 channel. It is 
reasonable to fix the wraparound contribution C. 
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TABLE VII: Summaries of the lattice results for the fitted energies of the nn system for the I — channel. The second 
block shows the energies in lattice units, where the first uncertainties are statistical, the second ones are the estimates of the 
systematic uncertainties. Column three shows the fitting range, and Column four shows the number of degrees of freedom (dof) 
for the fit. The six block shows the center-of-mass scattering momentum k 2 in GeV, and Column seven gives the product of 
pion mass and scattering length: m^o^°, where the first uncertainty is statistical and, the second one is systematic. 



Ensemble 


aE 


Range 


X 2 /dof 


A; 2 [GeV 2 ] 




(0.00484, 0.0484) 

(0.005,0.05) 

(0.007,0.05) 


0.33226(63)(78) 

0.3013(16)(18) 

0.3499(23)(26) 


9-24 
11 - 24 
9-32 


13.3/12 
14.6/10 
28.3/20 


-0.00572(21)(24) 
-0.00791 (71) (76) 
-0.0140(12) (13) 


0.476(25) (29) 

0.811(123)(133) 

1.181(202)(223) 
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FIG. 9: (color online). The "effective energy" plots as the 
functions of D m i n for the 7r7r scattering in the 7 = channel 
in lattice units. 



Nevertheless, one thing greatly comforting us is that 
the interaction among two pions in this channel is not 



too weak such that the discrete energies in a torus are 
shifted relatively bigger than that of the 1 — 2 channel 
from the values relevant for noninteracting pions, and 
as we can see from Tables [TT| and IVII1 the energy shift 
between the interacting and non-interacting 7T7t states is 
not a too small fraction of total energy. This indicates 
that the rigorous calculation of disconnected diagrams is 
at present the most important thing. 

The center-of-mass scattering momentum k 2 is calcu- 
lated by Eq. ([5]) with pion masses listed in Table UH and 
then the corresponding s-wave scattering lengths m w a^° 
can be obtained through Eq. ©. All of these values are 
summarized in Table IVII1 where the statistical errors of 
k 2 are calculated from the statistical errors of the ener- 
gies aE and pion mass am,, and its systematic errors are 
only estimated from the systematic errors of aE. Like- 



wise, the statistical for 



,1=0 



are computed from the 



statistical errors of k 2 and am^, and its systematic errors 
are estimated from the systematic errors of k 2 and one 
finite volume effect (27j . 

As already explained in the previous section, the de- 
pendence on the effective range r is small, and the range 
truncation actually leads to the finite volume correc- 
tion at 0(L~ e ) in Liischer formula (JB]) _26J. In prac- 
tice, we compute this correction for each lattice ensemble 
as suggested in Ref. [26||. This finite volume corrections 
have been combined in quadrature to the systematical 
errors listed in Table I VIII Other sources of systematic 
uncertainty like: nonuniversal exponentially suppressed 
corrections [HI, pion mass correction [H, Hl|, H^] , etc. are 
believed to be very small as compared with the rather 
large systematic error of the energies aE or we currently 
do not have enough computational resources to fulfil it. 
With more reliable calculation of the energies aE of the 
7T7T system in the I = channel in the future, these effects 
should be eventually incorporated into the more sophis- 
ticated lattice computation. 

We should point out that, in this work, we do not 
quote our results for the (0.01,0.05), (0.097,0.0484) and 
(0.0194, 0.0484) ensembles due to two considerations: 
First, the vacuum contributions of these ensembles are 
noisy (see Fig. [3]), and it is pretty hard to see the clear 
plateau (see Fig. ^ in the "effective energy" plots. Sec- 
ond, the presence of the a resonance is clearly presented 
in low energy (53l . |84| , and thus it should be necessary for 
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us to map out "avoided level crossings" between a reso- 
nances and 7T7T states with isospin-0 to secure the reliable 
scatterin g le ngth as investigated in the irK scattering in 
Refs. Luckily, as studied in Refs. [110111, the 

contaminations from a meson for three lattice ensembles 
with small pion masses are negligible. Therefore, we only 
consider these results in the rest of the analysis. 



VI. CHIRAL EXTRAPOLATIONS 



We should remark that the lattice calculations re- 
ported here used two lattice spacing of 0.15 fm and 
0.12 fm. Thus, it is meaningless to directly compare the 
energies aE of the tttt system. However, on the assump- 
tion that the Liischer technique properly explains for the 
finite volume dependence of the energies aE for these lat- 
tice ensembles, we can compare m^a^ 2 and m T a^° for 



two lattice spacings [27] , and we observe such agreement 
with statistical error in Tabic [V] 



In this work, we employed the rather small pion masses 
ranging from 240 MeV to 463 MeV, which are still larger 
than the physical one. Therefore, xPT is needed to carry 
out a chiral extrapolation of the scattering lengths to the 
physical pion mass. The resulting NLO xPT formulas, 
which can be directly built from the results in Ref. @ 
(see Appendix lAl for details), are described as [27l. |33| 



1=0 



ra^a 



1=2 



167T/2 



1 



91n^ 

J 7T 



1 + 



5 Kit (M — /ir,phy) 



16^ 2 /2 
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phy J 



(20) 



(21) 



where the values of and f n listed in Table |TT] are 
inserted into, and the X-PT renormalization scale is fixed 
at the physical pion decay constant /i = /^phy Where 
and whereafter a quantity with a "phys" subscript are 
referred to as the value of that quantity in the physical 
case. The l^°(fj,) and ^= 2 (/x) are the combinations of the 
LEC's in xPT at a quark-mass independent scale \i pBT - 
[27J • From the discussions in Appendix [5] the l^z {p) 
and l^ 2 (fi) arc connected to the LEC's l n as d HEf 
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It should be noted that the Eqs. (|20p and (f2"Tj) are ex- 
pressed in terms of the full /„■ computed on the lattice, 
and not the physical value /^.phy In reality, in the chiral 
expansion, the difference between utilizing /„. and f VlP hy 
in the argument of the logarithm only alters scattering 
lengths at higher orders [H |26|. 

As recommended in Refs. [25U271 ] . we will carry out 
the extrapolation of the products m^a^ 2 and m^a^ 
by means of the ratio / '/ w in place of m n . From Ap- 
pendix lAl we note that extrapolating in m„/ / w in lieu of 
ni-n does transform the representations for m n a^ 2 and 
m n a^° but only at NNLO or higher. Additionally, since 
m w // ff is a dimcnsionless quantity, there is no systematic 
error arising from the scale setting [25M27j . 



A. 1 = 2 channel 

We are now in a position to fit lattice results of m^a 1 ^ 2 
in Table |V] to the NLO xPT functional form {H} to ob- 
tain low energy constant l^ 2 (fi = / w , p hy), then the ex- 
trapolated value at the physical point (m w a^ 2 ) p hy S can 
be obtained. The lattice-calculated values of m^a]^ 2 as 
the function of m,/ f n are shown in Fig. 1101 and the outer 
error on the extrapolated result represents the systematic 
error and statistical error combined in quadrature. The 
one-loop xPT fit curve is displayed by the black solid 
line, and the red plus point indicates its physical s-wave 
scattering length: (m 7r a^ 2 )phys, which is the chiral ex- 
trapolation of the m^a^ 2 at the physical limit. In the 
same figure, we present the tree-level prediction as well. 
It is important to note that lattice data manifests pretty 
small displacement from the tree- level forecast. Addi- 
tionally, we notice that our lattice results for ro,B^ 2 are 
in general agreement with the one-loop formula. In fact, 
the deviation of (m w a^ 2 ) p hy S from tree-level prediction 
is a natural aftermath of NLO x?T fitting [27j |. 



I I 1 " 

— One loop 
-- Tree Level 
O Lattice data 




FIG. 10: (color online). The lattice-measured values of 
a function of m^j 'f n - The red plus point indicates 
the scattering length at the physical limit: (m^a^ 2 ) p hy S . 
The shaded bands correspond to statistical (inner-cyan) errors 
and statistical and systematic errors combined in quadrature 
(outer-yellow). The solid (black) curve is the central value of 
the NLO xPT fit. The dashed (magenta) line is the tree- level 
XPT prediction. 



In principle, we can fit our lattice-calculated data to 
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TABLE VIII: A compilation of the various theoretical (or phenomenological) , experimental and lattice QCD determinations 
of m n ai^ 2 extracted from the literature. Together with every reference, for an easier comparison, the first author name or the 
collaborations are given. The first uncertainty is statistical and second one is systematic if provided. 



Rcf. 


7=2 


Remarks 


This work 


-0.04430 ± 0.00025 ± 0.00040 


The calcualtion made in this paper. 


Yagi(2011) [31] 


—0.04410 ± 0.00069 ± 0.00018 


Extrapolation with NNLO x?T- 


Xu(2010) [27] 


-0.04385 ± 0.00028 ± 0.00038 


Using two flavors of maximally twisted mass fermions. 


NPLQCD(2007) [26] 


-0.04330 ± 0.00042 


Error combines statistical & systematic errors in quadrature. 


CLQCD(2007) [20] 


-0.0399 ± 0.0070 


The result from Scheme I of anisotropic lattices. 


NPLQCD(2005) [25] 


-0.0426 ± 0.0006 ± 0.0003 


With fully-dynamical domain- wall valence-quark propagators. 


Du(2004) [23J 


-0.0467 ± 0.0045 


Using anisotropic lattices in an asymmetric box. 


CP-PACS(2004) [24] 


-0.0413 ± 0.0029 


Compensating the mass dependence of the scattering length. 


JLQCD(2002) [22] 


-0.0410 ± 0.0069 


Selecting the result from EXP which employs a single exponential. 


Albaladejo(2012) [101 


-0.0424 ± 0.0012 


Employing unitary chiral perturbation theory. 


Guo(2009) [9] 


-0.0444 ± 0.0011 


Providing full results for all the contributing C(p 6 ) couplings. 


Sasaki(2008) [8] 


-0.0431 ± 0.0015 


Obtaining directly from the tvtt wave function. 


MILC(2004) [68] 


-0.0433 ± 0.0009 


Using MILC's determinations of LECs along with Roy equations. 


Zhou(2004) [7] 


-0.0440 ± 0.0011 


Chiral unitarization with crossing symmetry & phase shift data. 


CGL(2001) [3J 


-0.0444 ± 0.0010 


Two loop accuracy. 


Weinberg(1966) [1] 


-0.04557 ± 0.00014 


Tree level prediction. 


NA48/2(2011) [15, 16} 


-0.0429 ± 0.0044 ± 0.0016 


With independent experimental errors & different theoretical inputs. 


E865(2003) [11] 


-0.0454 ± 0.0031 ± 0.0010 


With the xPT constraints in the analysis. 



the NNLO x PT fo rm for m^a 1 ^ 2 @, H (see the concrete 
form in Eq. (JAToJ) ) as it is done in Ref. [3l|, since we 
have six lattice data at our disposal. In the meantime, 
we can make an estimate of NNLO LECs with the care- 
ful analysis for the chiral extrapolation of m^a 1 ^ 2 , since 
we have lattice data points at the lighter quark masses. 
However, the NNLO fit has larger errors in both l 1 ^ 2 and 
m-n-a^ 2 than those with NLO fit as shown in Rcfs. (25r - 

[27|,[3l|, and the errors of the LECs l [ ^ J=2 and if] J=2 
are rather large like the corresponding obtained values 
in Ref. [Hj]. Therefore, the calculations of NNLO LECs 
with physical meaning can not be obtained in this work. 
A rigorous NNLO xPT fit should wait for more lattice 
data at pion masses further closer to the physical point 
than we presently have. Admittedly, the resulting NNLO 
extrapolated value of m^aj'z 2 is indeed in harmony with 
NLO fit as we expect [27[. Actually, we use the NNLO 
xPT functional form to estimate systematic errors due 
to truncating the xPT series to NLO form (25| . 

As practiced in Ref. [27|], we only consider three ma- 
jor sources of systematic uncertainty on the extrapolated 
value of m^a^: 2 and l 1 ^ 2 - First, the lattice-calculated 
systematic errors of m^a^z 2 per ensemble are spread 
by the chiral extrapolation [271 ] . Second, the system- 
atic error inherently stems from NLO x?T fit (25), l26j |. 
which can be roughly calculated by taking the discrep- 
ancy between the NLO x?T extrapolated value from all 
six data sets and that from "pruning" the heaviest data 
set [27| ■ Third, the experimental errors on and / w (8(| 
give rise to another important source of systematic er- 



ror [271 ] . All three components are combined in quadra- 
ture to make up the entire systematic error. Taking the 
latest PDG data [86[ for the most accurate charged pion 
mass = m„+ = 139.57018(35) MeV and pion decay 
constant f n = = 130.41(20) MeV, where a couple 
of experimental errors are added in quadrature, hence 
"^tt/Ztt = 1.07024(166), we finally secure the upshots 

m^a 1 ^ 2 = -0.04430(25) (40); 
C= 2 (M = /., P hys) = 3.27(.77)(1.12), (24) 

where the first uncertainty is statistical and the second 
one is an estimate of the systematic error. 

These outcomes are comparable with the aforemen- 
tioned results of theoretical (or phenomenological) stud- 
ies [lill, HE! Gil; experimental determinations [HI, HE El] 
and lattice calculations [25T - I271 . l3ll ] within statistical er- 
rors. The relevant results for m T a^ 2 are courteously 
compiled in Table IVIIII The first group is lattice QCD 
results. The second one is theoretical (or phenomenolog- 
ical) studies. Also included are two experimental values 
in the third group. 

To make our demonstrations of these results more in- 
tuitive, they are offered graphically in Fig. [TTJ where we 
clearly note that the various results for every m^a 1 ^ 2 are 
fairly compatible with each other within errors. 

Our calculation of the LEC l 1 ^ 2 is satisfactory as well, 
although it is just about 25% precision, while it can 
be comparable with relevant results obtained by phe- 
nomenological predictions experimental determina- 
tions [HIIII, [lH, and lattice calculations [H-H^, [Mj]. 
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FIG. 11: (color online). A collection of various lattice QCD, 
theoretical (or phenomenological) and experimental results of 
m^a^ 2 listed in Table lylHl The lattice studies are shown by 
red diamonds, purple circles are theoretical (or phenomeno- 
logical) predictions, and the experimental ones are given by 
green crosses. Our result is indicated by a black square. For 
an easier comparison, the (cyan) inner strip corresponds to 
the statistical error whereas the (yellow) outer strip represents 
the statistical and systematic errors added in quadrature. 



The relevant values of l 1 ^ 2 are collected in Table IIXI 
The first group is lattice results. The second one is phe- 
nomenological and experimental determinations, which 
are transformed directly from the experimental and phe- 
nomenological results of m^a^ into ll=f at NLO xPT 
as preformed in Ref. [I?}- 



The reason why we made a significant improvement in 
precision over our previous work [33j is the recent com- 
prehension of various lattice-spacing artifacts (in partic- 
ular wraparound effect). In fact, approximate 0.5% accu- 
racy of our ultimate result for m^a^ 2 is typically joint 
efforts from lattice QCD and xPT- This can be under- 
stood from two aspects: First, we have lattice data closer 
to the physical point, which have relatively smaller un- 
certainties for m^a^ 2 '. Second, the chiral extrapolation 
of m^a 1 ^ 2 is considerably restricted by X-PT and m^a 1 ^ 2 
is solely predicted in terms of m T // T at LO and depends 
uniquely upon a LEC, l 1 ^ 2 , at NLO. This means that the 
statistical error of NLO xPT extrapolation of m^a^ 2 
solely rests on the statistical error of l^ 2 . Consequently, 
although our lattice-calculated results of m^a 1 ^ 2 are only 
with about 6% ~ 10% accuracy, we still obtain less than 
1% precise determination of {tyi- k q<- K tz )phys- 



TABLE IX: Some values of l^ 2 extracted from the literature. 
The first uncertainty is statistical, and the second one is sys- 
tematic if present. The first group is lattice QCD results. The 
second one is phenomenological and experimental determina- 
tions, which are transformed directly from the corresponding 
results of m^a 1 ^ 2 into l 1 ^ 2 at NLO xPT [13] . 



Ref. 


l'= 2 

7T7T 


This work 


3.27 ± 0.77 ± 1.12 


Yagi(2011) [31] 


5.8 ±1.2 


Xu(2010) [27] 


4.65 ± 0.85 ± 1.07 


NPLQCD(2007) [26] 


6.2 ± 1.2 


NPLQCD(2005) [25] 


3.3 ±0.6 ±0.3 


CGL(2001) [3J 


3.0 ±3.1 a 


NA48/2(2011) [15,26] 


7.5 ±13.3 ±4.8 6 


E865(2003) 21] 


0.0 ±9.4 ±3.0 



a It is interesting to note that if we make use of Eq. J23D with 
the values of ij reported in Ref. 0], and necessary PDG values, we 
obtain i£= 2 = 2.0 ±3.1. 

' Using the data from Ref. [l5l | . 



B. 



channel 



We are now in a position for fitting lattice results of 
m n ai=° in Table EH to NLO xPT functional form j20]) 
to secure the low energy constant l^°(fi = / w , p h y ), 
then obtain the extrapolated value at the physical point 
(m 7r a^°)phy S ' The lattice- measured values of m ff a^° 
as the function of m„/ 'f„ are demonstrated in Fig. [T2l 
The one-loop x?T fit curve is displayed by the black 
solid line, and the red circle point indicates its physical 
s-wave scattering length: (TO 7r a^ ) p hy S , which is the chi- 
ral extrapolation of the m^a^ at the physical point. 
In the same figure, we present the tree-level prediction 
as well. It is important to note that lattice data mani- 
fests a rather large displacement from the tree-level xPT 
forecast, which is well consistent with one conclusion 
of Ref. [|| that the NLO corrections make about 25% 
modification to tree-level prediction. Additionally, we 
can notice that our lattice results for m^a^ generally 
agree with the one-loop formula. The large deviation of 
(m 7r a^°) p hys from tree-level prediction is totally a nat- 
ural aftermath of the NLO xPT fitting. 

In this work we can not fit our lattice-calculated data 
to the NNLOxPT functional form ([20]) [3, l3j (see con- 
crete form (|A10[) ). since we just have three lattice data 
at our disposal. A NNLO xPT determination should 
wait for more lattice data at pion masses further closer 
to the physical point than we now have. Admittedly, 
as explained in detail in Ref. [||, the NLO correction 
increases the LO prediction by about 25%, and NNLO 
correction further raises the LO prediction more than 
10%, this mean that the NNLO xPT determination of 
should be significantly more rigorous than the 



ra^a 



1=0 



NLO X-PT determination of (m w a^ 7r ) p hys, and this can 
in part explain the relative large error for our final NLO 
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TABLE X: A compilation of the various theoretical (or phenomenological) , experimental and lattice QCD determinations 
of m^-a^ extracted from the literature. Together with every reference, for an easier comparison, the first author or the 
collaborations are given. The first uncertainty is statistical and the second one is systematic if given. 



Ref. 


ATI 7f Ch'jY 7f 


Remark 


This work 


0.214 ± 0.004 ± 0.007 


Full QCD. 


Fu(2012) [331 


0.186 ± 0.002 


Partially quenched QCD. 


Albaladeio(2012) [101 


0.219 ± 0.005 


Employing unitary chiral perturbation theory. 


Guo(2009) [9] 


0.220 ±0.005 


Full results for all the contributing 0(p b ) couplings. 


Caprini(2008) [88] 


0.218 ±0.014 


Using a large class of analytic parameterizations. 


Yndurain(2007) [87] 


0.233 ±0.013 


Extrapolating to the pole of the sigma resonance. 


Zhou(2004) [7] 


0.211 ±0.011 


Using chiral unitarization with crossing symmetry & phase shift data. 


CGL(2001) [3] 


0.220 ±0.005 


Two loop accuracy. 


Weinberg(1966) [1] 


0.1595 ±0.0005 


Tree level prediction. 


NA48/2(2011) [15, i6] 


0.2210 ± 0.0047 ±0.0015 


With independent experimental errors & different theoretical inputs. 


E865(2003) [11] 


0.216 ±0.013 ±0.002 


With the xPT constraints in the analysis. 




m /f 

FIG. 12: (color online). The lattice-measured values of 
m n a^° as a function of m n /f n . The red circle point indicates 
the scattering length at the physical limit: (m^a^ ) p hy a . 
The shaded bands correspond to statistical (inner-cyan) errors 
and statistical and systematic errors combined in quadrature 
(outer-yellow). The solid (black) curve is the central value of 
the NLO xPT fit. The dashed (magenta) line is the tree- level 
XPT prediction. 



%PT extrapolated result of (m 7r a^ ) p i 1 y S - As a conse- 
quence, the systematic error from truncating the xPT 
series to the NLO form should be considered. 

In this work, wc only consider two major sources of 
systematic errors in the extrapolated value of m,o^ 
since the systematic error from the experimental errors 
on TO^ and f n is found to be pretty small as compared 
with its statistical error. First, the lattice-computed sys- 
tematic errors of m^a^ per ensemble are propagated 
via the chiral extrapolation. Second, the systematic er- 
ror inherently stems from NLO fit, which is computed 
by taking the difference between the extrapolated values 
from NLO fit to all three data set and that from 'crop- 
ping" the heaviest data set [27} . All two parts are added 



in quadrature to give the whole computed systematic er- 
ror. Wc get the final upshot 

mW=* = 0.214(4)(7); 
l^°(n = f ntphys ) = 43.2 ±3.5 ±5.6, (25) 

where the numbers in the first and second parentheses are 
the statistical and systematic uncertainty, respectively. 

These results can be fairly comparable with the above- 
mentioned results by theoretical (or phenomenological) 
studies 0, HES HIH (except the tree-levehpredic- 
tion), and experimental determinations [ill flH ITq]. The 
relevant results for m T o[j are compiled in Table IXl The 
first group is lattice results. The second one is theoret- 
ical (or phenomenological) studies. Also contained are 
two experimental values in third group. 

To make our report of these results more intuitive, 
these results are given graphically in Fig. [T3] as well, 
where the various results of m w a^° are compatible with 
each other within errors except the tree-level predic- 
tion [l| and our crude lattice study [33j] . 

Our calculation of the LEC is satisfactory as well, 
although it is about 5% accuracy, nevertheless, as we 
show soon, this result can be comparable with the re- 
sults of phe nomenological 0] and experimental determi- 
nations [II [IS El, and lattice studies HSH^iJ. The 
relevant values of are collected in Table IXll The first 
group is lattice results. The second one is phenomeno- 
logical and experimental determinations, which are con- 
verted directly from the experimental phenomenological 
results of m T a^° into at NLO x?T as conducted 
by Xu et al. in Ref. [13 for the J = 2 channel. 

The remarkable improvement in accuracy over our pre- 
vious results [33j] is an joint effort from lattice QCD and 
XPT First, we have measured the lattice data closer to 
the physical point, which have smaller uncertainties. Sec- 
ond, chiral perturbation theory considerably constrains 
the chiral extrapolation of the product m T a^°, which is 
uniquely predicted in terms of ra^j f n at LO and relies 
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FIG. 13: (color online). A collection of various lattice QCD, 
theoretical (or phenomenological) and experimental results of 
m n a^° listed in Table IXl The red diamonds are lattice deter- 
minations, purple circles are theoretical (or phenomenologi- 
cal) studies, and experimental ones are represented by green 
crosses. Our result is shown by a black square. For an easier 
comparison, the (cyan) inner strip corresponds to the sta- 
tistical error whereas the (yellow) outer strip represents the 
statistical and systematic errors added in quadrature. 



TABLE XI: Some values of extracted from the litera- 
ture. The first uncertainty is statistical, and the second one 
is systematic if present. The first group is lattice simulation 
results. The second one is phenomenological and experimen- 
tal determinations, which are directly transformed from the 
corresponding results of m^a 1 ^ into at NLO x?T 



t i=o 

L 7Y7Y 


This work 


43.2 ± 3.5 ± 5.6 


Fu(2012) [33] 


18.7 ± 1.2 


CGL(2001) [3] 


48.5 ±4.3 a 


NA48/2(2011) [15] 


49.3 ±4.1 ±1.3 


E865(2003) [11] 


45.0 ± 11.2 ±3.5 



a If we make use of Eq. (122[ i with the values of ij reported in 
Ref. 0, and required PDG values, we get i£=° = 32.4 ± 2.3. 



solely on a LEC, l„^°, at NLO. This suggests that the 
statistical error of NLO xPT extrapolation of m ff a^° 
solely rests on the statistical error of Consequently, 
although our lattice-calculated results of m^a^ are only 
with 5% ~ 17% accuracy, we still obtain about 2% pre- 
cise determination of (m 7r a^ ) p hys- 



VII. SUMMARY AND OUTLOOK 

We have reported results of a lattice QCD calculation 
of the s-wave 7T7t scattering lengths for both 1 = and 
2 channels on the MILC "medium coarse" (a w 0.15 
fm) and "coarse" (a ~ 0.12 fm) lattice ensembles with 



the 2+1 flavors of the Asqtad-improved staggered sea 
quarks. We exploited the moving wall sources without 
the gauge fixing 1^ 19l| t o compute all the four diagrams 
assorted in Refs. jTa Tl8| . and viewed a clear attractive 
signal for the 1 = channel and a good repulsive one for 
the 1 = 2 channel, respectively. Moveover, extrapolat- 
ing our lattice data of the s-wave scattering lengths for 
both isospin eigenstates to the physical pion mass gives 
the scattering lengths: m^a 1 ^ 2 = -0.04430(25) (40) and 
771^(4=° = 0.214(4)(7) for the I = 2 and channels, re- 
spectively, which are in fair agreement with the current 
theoretical (or phenomenological) predictions to one-loop 
levels and the present experimental reports, and can be 
comparable with other lattice studies. 

After our extremely crude estimation of the 7T7t scat- 
tering length in the 1 = channel in Ref. J33|, and this 
relatively more sophisticated computation, we can fairly 
claim that even with the limited computing resources, 
the lattice calculation of the 7T7t scattering length in the 
1 = channel is feasible although this work is absolutely 
needed to be further improved, and the various sources 
of systematic error yet need to be clarified thoroughly. 
Most of all, from this work, we found that the rule of 
thumb estimation of lattice ensemble with the Goldstone 
pion mass eligible to study the s-wave 7T7t scattering in 
the 1 = channel should be less than about 300 MeV 
(the smaller, the better), which is very helpful for the 
people to pursue this fascinating enterprise. We view it 
as one of the important results of this work. 

As we revealed, a reasonable signal can be gained for 
the (0.00484,0.0484) and (0.005,0.05) ensembles in the 
vacuum diagram of the 7T7r scattering. In principle, the 
signal-to-noise ratio can be further enhanced by launch- 
ing the same calculation on the lattice ensembles with 
a smaller pion mass (of course we can also improve the 
statistics by using more lattice gauge configurations or 
performing the calculation on a larger volume). In ad- 
dition, the behavior close the physical point is intensely 
influenced by the chiral logarithm term, so an extraction 
of the 7T7r scattering lengths without a long extrapolation 
is still highly needed to guarantee the convergence of the 
chiral expansion. Fortunately, the MILC collaboration 
has generated enough lattice ensembles whose Goldstone 
pion masses are smaller than or close to 240 MeV [III [9(3] 
(e.g., the fine (0.00155,0.0310) ensemble, whose Gold- 
stone pion mass is about 177 MeV). Furthermore, as we 
explained early, the NNLO x?T will be the proper phys- 
ical functional form to fit the lattice data (at least four 
data point) of the I = itit scattering length and it needs 
more lattice data near the physical point. I have an im- 
petus to do these works. However, it is beyond the scope 
of this work since this will need an astronomical amount 
of computing allocations. We will enthusiastically appeal 
for all the possible computational resources to accomplish 
this peculiar and challenging task. 

It is well-known that 7T7t scattering in the 1 = channel 
is challenging and stimulating phenomenologically due to 
the existence of the a resonance. In the work, we have ex- 
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hibitcd that the tttt scattering for the I = channel can 
be reliably calculated by the moving wall sources without 
gauge fixing [IH, HH . It gives us an anticipation that this 
technique can be successfully exploited to study the a 
resonance. In our previous work [771 ]. we have evaluated 
the a mass from lattice QCD, and found that the decay 
a — > tttt is allowed kincmatically only for enough small u 
quark mass. This work and our lattice calculation for the 
tttt scattering lengths delivered in this paper will encour- 
age the researchers to study the a resonance. We have 
been investigating the a resonance parameters with the 
isospin representation of (I,I Z ) = (0,0), and the prelim- 
inary lattice results are already reported in Ref. [5J] . 

Additionally, for tttt scattering in the 7 = channel, 
we realize an important issue that the presence of a res- 
onance is possible in the low-energy, and thus it should 
be necessary for us to employ the variational method [HI 
to secure the rigorous scattering length as investi gate d in 
the ttK scattering in the 1 = 1/2 channel in Refs. [53ll84| . 
Since we only make use of the relative small quark masses 
to study the tttt scattering in the 1 = channel, we can 
temporarily and reasonably overlook this contamination 
in the present study as already explained in Ref. [53| . 
However, we should bear in mind that this issue should 
be settled in the more sophisticated lattice examination. 
It will be very interesting to systematically study this 
pollution for the 7T7t scattering in the 1 = channel. 

Admittedly, due to the intense theoretical and exper- 
imental efforts put into the scalar-isoscalar and scalar- 
isovector sector of the meson-meson scattering recently, 
studying the KK scattering length on the lattice is grad- 
ually beco ming a very interesting enterprise. As pointed 
out in Ref. [9l[ , the robust calculation of the tttt scatter- 
ing lengths (in particular for the 1 = channel) will nat- 
urally encourage us to study other challenging systems 
like KK, etc. Physically, as explained in Ref. pQ, l65l. l92| . 
studying KK is very important, and the calculation of 
the s-wave KK scattering length in the 1 = channel 
is a genuine two coupled-channel problem [65l [93j , where 
the system can be approximately described by only tttt 
and KK two channels (we refer to tttt as channel 1 and 
KK channel 2), then the 5-matrix is a 2 x 2 unitary 
matrix which contains 3 real parameters (6f| [93| (the s- 
wave KK scattering in the 1 = 1 channel can be treated 
analogously HH , and the lattice study of ttij scattering is 
then highly desirable.). Therefore, it is absolutely nec- 
essary to incorporate the s-wave I = 0, 7T7r-channel for a 
physical calculation of the s-wave KK scattering length 
in the 1 = channel. The generalized Liischer's formula 
in this case gives a relation among these three parame- 
ters, all of which are functions of the energy [93J. since 
some of these parameters are still poorly measured in the 
present experiments, the lattice calculation is valuable 
and highly desirable. With our lattice efforts on channel 
2 in Ref. [6(|, at present, if we can compute the s-wave 
tttt — > KK scattering in the 1 = channel, in principle, 
we can solve this problem. We are launching a series of 
lattice studies for these aims. 
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Appendix A: Scattering length of rnr in xPT a ^ 
NNLO 



In Ref. [331 ] , we provided the compact continuum three- 
flavor xPT form for the s-wave tttt scattering length for 
isospin-0 at the NLO in the continuum limit of QCD by 
constructing from Appendix C in Ref. 0] . Here we follow 
the original derivations and notations in Refs. 0-QL Hl[ 
to derive its compact form at the NNLO. 

The tttt scattering lengths are provided at the NNLO of 
xPT in Refs. 0-Q|, and the s-wave tttt scattering length 
in the 1 = channel is described as 0-QI 



ra^a 



1=0 



7mt 



32tt/2 I " ' 7 



1 + - [49 + 5&i + 126 2 + 486 3 + 326 4 ] 
7045 215tt 2 



63 
192- 

-64& 4 + — h 



126 



106i + 24& 2 + 966 3 



+ 0(x 4 ), 



(Al) 



10 The numerical calculations of this work are unceasingly carried 
out for more than two years. We should especially thank Prof. 
Hou qing's continuous encouragements and comprehensive sup- 
port, without his kind and selfless help, it is not possible for us 
to launch this work, and have an opportunity to do it. 
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where x = m 2 /(167r 2 / 2 ) is the chiral expansion param- 
eter and 6j's are dimensionless combinations of the cou- 
pling constants introduced in Refs. 0,0] to parameterize 
the pion scattering amplitude. After some strenuous al- 
gebraic manipulations, we can check that 



5&i + 126 2 + 486 3 



32 6 4 = b ll-™xL 2 
2 6 
20- 40- 5~ 1 - 63 

~ ( 388 r 472 r orr ir - 1405 \ 
+xL f — — k - — Z 2 - 35Z 3 + 154Z 4 + — J 

+x JjW 4 + ^[ 2 f 4 - 15^4 + 35Z~ 4 2 - ^~ 3 2 

364 r 1336 r 141- , ~ 162719 

+ ~ h + — h + ~ h - 126/4 + ^32- 
373 
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5?~i + \2r 2 + 48f 3 + 32fz 



(A2) 



where the low-energy constants li, fi are the quark mass 
independent couplings from the subleadingorders £4, 
of the effective Lagrangian, respectively [3j, and renor- 
malized at the running scale /1, and L = ln(/x 2 /m 2 ). 

Inserting Eq. (|A2|) into Eq. (|A1[) and considering ex- 
pression B. 3 of Ref. [|[ 



b 85 f 2 
65 - 72 L 

7 



r ,'7r 107 r 6251 
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we recast the result in the order of x as 
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The right-hand side (rhs) of Eq. (|A1[) is scale indepen- 
dent |3|. On the whole, the rhs of Eq. (|A3|) is scale in- 
variant as well. Therefore, in principle, we can select the 
running scale [i stochastically. However, when fitting our 
lattice-obtained scattering lengths as a function of x, it 
is highly desired for us to settle the fitting parameters 
quark-mass independent. As it is done in Ref. [3ll ]. we 
select /1 = 47r/ 7 r ! ph y . Using this scale, we can see 



L(/j, = 47r/ 7r , p hy) 



lnx - 2xk(^i = Airf^^y) 



+2x\nx + 0(x 2 ), 



(A5) 



where we exploit the chiral expansion of the pion decay 
constant = /^phyjl + xk + 0(x 2 )} 0. 

Plugging Eq. (|A5|) into Eq. (|A3[) . and rearranging the 
result in the order of x, we achieve the tttt scattering 
length in the 1 = channel as 
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The continuum xPT forms for the s-wave tttt scatter- 
ing length in the 1 = 2 channel aj^ 2 at the NNLO was 
presented by Yagi et. al. in Ref. [3l| . namely, 
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In the above equations, is the pion decay constant, 
which is originally written as (around 92.4 MeV) 0- 
4]. In the present work, a/2-FV is denoted by /„■ (about 
130 MeV) for the convenience of the fitting our lattice 
data. Then the above equations are recast as 
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where l^ls are the combinations of LEC's in xPT at a 
quark-mass independent running scale since all the LEC's 
are independent of quark mass, therefore, we can regard 
them as the fitting parameters in the chiral extrapolation 
of the s-wave -kit scattering lengths [3l[ . 

From Eqs. (|Ai]t . (jA6|) . (|A7]) and (|A8j) . we can easily 
get its specific forms of and Z^ 2 , which are related 
to the Gasser-Leutwyler coefficients li as Q, HH 
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where we consider the equality /„ = l n + ^ n ( m n/l j2 ) 0- 



These are the forms what we used in our previous 
work (33|. For the other l^ls, its explicit forms are given 
or can be inferred from Eqs. (|A4|) , ()A"6|) , (|AT| and (|A8[) . 
It is interesting and important to note that if we opt the 
running scale fi = /^phy, we obtain the same expressions. 
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